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I In this work, we study the hght stop pair signals at the LHC. We explore the 

O ■ SUSY parameter space with non-universal gaugino and third generation masses at 
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Q . the GUT scale. Recent LHC SUSY search results based on 35pb~'^ and lfb~^ of data 

are implemented to put the limits on stop pair events. The dark matter relic density 
' and direct detection constraints are also taken into account. Detailed simulations 

on the signals and background for some benchmark points are performed, and it 
, is found that the stop pair signals usually escape the LHC search if the present 

_ cut conditions are used. We also explore the potential and sensitivity of ILC to 
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f"*) ■ probe such scenarios. It is found that the ILC can detect them with an integrated 
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I. INTRODUCTION 

Supersymmetry (SUSY) is one of the most attractive extensions beyond the standard 
model (SM) which offers a solution to the hierarchy problem and grand unification in gauge 
interactions. Moreover, the R-parity conserved SUSY models naturally provide a lightest 
supersymmetric particle (LSP) which is neutral and stable and can be a dark matter can- 
didate badly needed in order to interpret the cosmological observations. To guarantee the 
electroweak symmetry breaking naturally, the size of SUSY soft breaking terms is believed 
to be of the electroweak scale, of which the light Higgs boson mass is typically less than 
about 140 GeV and the masses of sparticles are around a few hundred GeV which are within 
the reach of LHC. The running of LHC can test supersymmetry by discovering a light Higgs 
boson and SUSY particles. 
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Search for supersymmetry is one of the prime targets of LHC Certainly the strong in- 
teraction sector of the SUSY models is the most important place to discover SUSY due to 
the large cross sections for the production of colored particles, such as squarks and gluinos, 
at the hadron colliders. As one of typically signature of such processes, the LSPs appearing 
in the final states lead to a large missing energy (experimentally a large missing transverse 
momentum). Therefore SUSY signature should show up best in the jets plus large missing 
energy channel. Currently, the most stringent bound is derived from this channel, which is 
the dominant signature for gluino and squark production. With lfb~^ of data at the LHC, 
no SUSY signatures are detected. The bound on the gluino and squark masses are set to 
be larger than about 800 GeV for the mSUGRA model [1, 2]. The recent release of LHC 
measurements motivated quite a few new works on SUSY [3] , where the implications of CMS 
and ATLAS results to various SUSY models are discussed. 

If nature chooses SUSY as its working principle, the null results in the supersymmetric 
searches at hadronic colliders (both Tevatron and LHC) may indicate that SUSY may hide 
from our probe in some ways. For example, one way is that the SUSY might be splitting 
in its color sector, i.e. color sparticles are very heavy (say more than a few TeV) and their 
cross sections are highly suppressed by masses. Then SUSY can only show up at LHC 
via processes from its electroweak sector. In order to find SUSY, we have to accumulate 
more data. In some worse cases, SUSY may even be missed at LHC. For example, almost 
degeneracy of the masses of chargino and neutralino leads to soft leptons which can escape 
the detection at LHC [4]. If it is such linear electron-positron collider or a muon 

collider should be necessitated in order to unravel such a scenario. 

Another possible way is that SUSY signature is buried deeply in the background events 
and can not been selected out with the current search approach. For example, in the 
compressed SUSY models [5, 6], where the mass difference between the squark and LSP is 
small, only few soft jets can show up in the final states. Thus the signature is so similar to 
the background events that the currently SUSY search selection conditions can not separate 
the signal events. 

The scenarios in which gluino and squarks of the third generation are light but degenerate 
with LSP may also be hidden to the LHC search [7-9]. One of such a possibility is the 
light stop scenario [10-27]. The stop can be the next-to-lightest supersymmetry particle 
(NLSP), which can be quite naturally realized in SUSY models, like mSUGRA. The large 
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top Yukawa coupling possesses two-fold effects to the stop quark mass spectra. First, the 
off-diagonal trilinear term, which is proportional to top Yukawa coupling, can lead to the 
largest mass splitting in squark sector and consequently produce a light stop as the NLSP. 
Second, since the lightest stop is mostly right-handed, the mass of the right-handed stop, 
m^^, is significantly reduced by the top Yukawa coupling via the renormalization group 
running. Light stop scenario is also well-motivated to explain the dark matter relic density 
measured by WMAP when the mass splitting between it and the dark matter candidate (the 
LSP) neutralino (x?) is small enough [13]. Such a case is also dubbed as the stop-neutralino 
coannihilation scenario. Furthermore, the electoweak baryogensis in the framework of MSSM 
favors a light stop [11, 12], i.e. the light stop can generate the first order phase transition 
which prevents the baryon asymmetry of the universe from being washed out. In order to 
guarantee this first order phase transition, the mass of the light stop should be light. 

In the mSUGRA scenario, all the gaugino and scalar masses are usually assumed universal 
at the GUT scale. However, such an assumption is not guaranteed by symmetry. If the F- 
term of the gauge kinetic function in the SU(5) SUSY GUT is not singlet, the gaugino 
masses at GUT would not be universal [28, 29]. The non-universal gaugino masses can also 
be predicted in the supersymmetric SU{5) x SU {?>) Hypercoior niodel proposed to solve Higgs 
doublet-triplet splitting problem in SU(5) GUT [30, 31], and in the supersymmetric partial 
unified model SU (4)c x SU{2) l x SU (2)/j[32]. For the scalar sector, notice the most stringent 
constraints on SUSY FCNC and CP violation processes only relate to first two generations 
sfermions, it is well motivated to consider the so-called "inverted scalar mass hierarchy" 
scenario [33-36]. In such scenario, the large masses of first two generations sfermions can 
solve the SUSY flavor and CP problems, while the third generation sfermions are still light 
to satisfy the naturalness conditions. These non-universal soft breaking parameters will 
change the running behavior of the RGE, and induce different sparticle mass spectra and 
search strategies from the mSUGRA. 

In this work, we will study a SUSY scenario with non-universal gaugino and third gen- 
eration masses at the GUT scale, and then explore the signatures of the light stop pair 
production (for some relevant studies, see [37-42]). We first scan the SUSY parameter space 
in the non-universal SUSY model. We have considered all the constraints on the SUSY 
parameter space. Present 1 fb~^ of LHC data at ATLAS and CMS are used to put limits on 
stop pair events. We also simulate the signatures in some benchmark points at the LHC in 
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detail. In general, the stop pair productions depend on stop mass parameter, while the light 
stop decay modes only rely on light slepton, chargino and neutralino mass spectra. Our re- 
sults can be easily extended to the scenarios with decoupled gluino and first two generation 
squarks. For the stop-neutralino coannihilation scenario to give correct dark matter relic 
density we find the SUSY signatures are hidden by the present cut conditions. We further 
present a study of the signals at the future linear collider machine. 

This paper is organized as follows. In Section 2 we start with the generic bounds on SUSY 
from LEP, Tevatron, and LHC, while we concentrate on the bounds to the light stop from 
its various decay channels searched by experiments. In Section 3 we analyze the dark matter 
bounds to the parameter space of non-universal SUSY models and study the mass spectrum. 
In Section 4 we discuss the recent LHC bounds for our light stop scenario, and select four 
bench mark points, demonstrate their mass pattern and decay features, and elaborate how 
the signature from these bench mark points can hide from the current search at LHC. We 
also study the sensitivity of ILC to these bench mark points. Section 5 is for discussions 
and conclusions. 

II. CURRENT EXPERIMENTAL BOUNDS ON LIGHT STOP SCENARIO 

In this section we describe the available bounds from the stop pair search at LEP, Tevatron 
and LHC in details. Several stop decay modes are investigated in [43], i.e. 1) t — ?■ xfb, 2) 
t — > uib, 3) t — > i'^ub, and 4) t — XiC- The recent bounds listed below depend on the mass 
spectrum and the decay mode of the sparticles. 

• ti — )■ bx~^ — > biuxi- The pair production of stop decaying via ti -> bx~^ has been 
investigated by the CDF collaboration with an integrated luminosity of 2.7 fb~^ [44]. 
The stop masses between 128 and 135 GeV are excluded at 95% independent of 
the branching ratio of x^ ~^ ^^Xv For m(xi) = 45GeV, 'm(xf) = 125.8GeV and 
Br{x^ iuxi) = 1, lower limit for stop mass can be set at 196 GeV. 

• ti — )■ biu: The most stringent bounds for stop pair production with this decay channel 
are given by the DO collaboration [45] with an integrated luminosity of 5.4 fb~^. 
The main final states are focused on bbe^fi^uu. The sneutrino is assumed to be 
the LSP or decay invisibly into uxi- The analysis is optimized for the mass region 
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Am = m^-^ — mc, = AOGeV or above. Stop masses up to 240 GeV are excluded for 
sneutrino masses around to 45 GeV, and sneutrino masses up to 120 GeV are excluded 
for stop masses around to 180 GeV. 

• ti — )■ tXi'- The recent search for the top partner T' (it can be the t' of the fourth 
generation model, the new heavy quark in the little Higgs model with T-parity, or 
the scalar top quark in the SUSY model) and T' — tx via the semi-leptonic mode 
pp —J- iuebqq'b + XX fuU-hadronic final states pp — qiq2hq3,<id> + XX by CDF 
collaboration are reported from the Ref. [46] and [47], with integrated luminosities of 
4.8 fh~^ and 5.7 fh~^ respectively. It is shown that the mass of T' can be bound up to 
360 GeV for mx < lOOGeV and 400 GeV for nix < 70GeV, respectively. When both 
top quarks decay semi-leptonically, the final states are the same as the first two decay 
modes, i.e. ti — ?■ bx^ and ti — )■ biu. The discovery of this decay channel with both the 
discovery hadronic and semi-leptonic channels can distinguish this decay mode from 
the others. 

• ti — cXi- If the the above processes are all kinetically forbidden, the loop induced 
flavor changing process ti — cXi might be dominant decay channel. The similar pro- 
cess ti uxi is always suppressed by the CKM matrix when ti — ?■ cXi kinematically 
opens. This channel is difficult to be detected if the mass splitting between ti and Xi 
is smaller which leads to two soft jets in the final states. The CDF collaboration had 
performed the research for this process with 2.6 fb~^ of data [48]. At least one of the 
jets is required to be tagged from a heavy-fiavor quark. The analysis is optimized for 
the mass region Am = m^-^ — mj> = AOGeV or above. Stop masses up to 180 GeV are 
excluded for neutralino masses around to 95 GeV. In addition, the results from LEP 
had excluded stop masses up to 90GeV for ti — cXi = 1. A more recent experimental 
search at DO collaboration can be found in [49]. 

• ti — > bff'xi'- For small mass splitting between stop and the LSP, the four body process 
ti — 7- bff'xi could be also important [50]. Typically, the objects in this channel are 
softer comparing with those from ti txi- 

• R-hadron : If the life-time of stop is long-enough due to the extremely small mass 
splitting or weak interaction between stop and the LSP [15], stop can form a bound 



6 



state R-hadron in the process of hadronization before its decay. If such R-hadron 
carrying electric charge, it might be observed in the inner tracker and even outer 
muon detector. Recently, CMS collaboration reported the limits for R-hadrons from 
pair production of stable stops based on 1.09 fb~^ of data [51]. Lower limit for stop 
mass can be set at 620GeV and 515GeV, corresponding to whether R-hadrons can leave 
observable signatures at the muon detector or not respectively. The search strategy 
of long-lived stops in MSSM is studied at LHC and can be found in the Monte Carlo 
study by Ref. [17]. 

III. ANALYSIS OF THE PARAMETER SPACE 

A. Theoretical scenarios 

In the mSUGRA scenario, all the gaugino masses are set to be universal as mi/2 at the 
GUT scale. However, this is a convenient assumption rather than a theoretical requirement. 
A non-universal gaugino mass sector is well motivated in many superstring or SUSY GUT 
models [28, 29, 37-39]. In the SUSY GUT, if there exists a holomorphic function /($) in 
the gauge kinetic function, gauginos will acquire masses via non-zero F-term of the /($) 

j d'9 f{<^),,W''W' + h.c. ~ ^(F(<I>)),,A'^A^ (1) 

where $ is a chiral field in the hidden sector related to SUSY breaking, {F{^))ab transforms 
in the Adj ® Adj representation of the underlying gauge group containing gauginos A"^, A*. 
For SU(5) SUSY GUT, {F{^))ab belongs to 24 O 24 = 1 © 24 © 75 © 200. Only if 
is a singlet as {F{^))ab ~ c6ab, gaugino masses are universal as in the mSUGRA. The 
non-universalities are generated when {F{^))ab belongs to other high representations. For 
example, the gaugino mass relationships at GUT scale are 2 : —3 : — 1, 1 : 3 : —5 and 
1 : 2 : 10 for /($) belonging to 24, 75, 200 respectively [28, 29]. Moreover, if different /($) 
representations appear simultaneously, other gaugino mass relations can be achieved. 

In the scalar sector, a general Kahler potential could also leads to mass non-universality 
[37, 38] 

J dVe KoQ^Q + KiQlQ^ + ... ~ Q,</.V, (2) 
where Kq, Kj, ... are real function of $. In the mSUGRA scenario, all the scalar 
fields acquire same mass under the universal Kahler potential assumption Kj ~ K6j. This 
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assumption is useful to suppress FCNC effects which are stringently constrained by experi- 
ments. However, solving the SUSY flavor and CP problems only require large masses of first 
two generations sfermions, while the third generation sfermions can still be light to satisfy 
the naturalness conditions [33-36]. 

Therefore, in this work, we treat these soft breaking mass parameters as free parameters 
at the GUT scale without imposing specific relations among them which can be derived from 
their underlying non-universal models. For the sake of simplicity, we choose the following 
seven input parameters as free parameters in our analysis: 

Ml/2, Ml/2,3, rrio, mo,3, ^o, tan (3, sign{fi) . (3) 

Compared with mSUGRA, only two extra input parameters at GUT scale, i.e. the gluino 
mass Ml/2,3 and the third generation sfermion mass mo,3, are added. 

B. Parameter Space Scan and Experimental Constraints 

In this subsection we scan the SUSY parameter space with the 7 free parameters at 
the GUT scale discussed above. The low energy spectra are calculated by solving the 
RGBs. In the scanning the GUT and the electroweak symmetry breaking scale are set to 
be Mgut = 2 X 10^^ GeV and Mewsb = y/^^HJ^) respectively. We consider the case 
where gluino is lighter than wino and bino at Mqut and put upon a constraint by requiring 
lOOGeV < Ml/2,3 < Ml/2 < SOOGeV. To obtain a lighter stop or stau which is necessary in 
order to yield a suitable DM relic density, we further assume that third generation sfermions 
are lighter than the other scalars in the first two generations, and allow them to vary in the 
range lOOGeV < mo,3 < niQ < 2000GeV. The trilinear coupling Aq and the ratio of vacuum 
expectation values tan /3 are chosen in the range of —1 < Ao/mo < 1 and 2 < tan/3 < 50, 
respectively. The sign of /i is taken to be positive, which is favored by several experimental 
constraints on 6 — > 57 etc. We utilize SuSpect [52] to calculate the SUSY particles spectra by 
solving the two-loop SUSY renormalization group equations. The top quark pole mass can 
affect the sparticle spectra and consequently modify dark matter relic density significantly. 
We take the top pole mass as mt(pole) = 173.1 GeV. 

Several phenomenology and astrophysics experimental constrains are taken into account 
in our scanning. Important flavor physics constraints: Br{b — )■ S7) = (3.55 ± 0.24) x 10~^ 
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[53], Br{Bs f-i+i^-) = (0 ± 1.4) x IQ-^ [54], 5r(5„ ^ tu)/SM = 1.28 ± 0.38 [53], are 
realized. As a conservative analysis, we only demand that the SUSY contributions pass 
these constraints to a 3a level. Another remarkably important constraint is from the muon 
anomalous magnetic moment — 2 measurement [55], which is adopted here as a bound 
—11.4 X 10"^° < — 2 < 9.4 x 10~^ as used in Ref. [40]. The lighter Higgs boson is required 
to be heavier than 114 GeV, while the mass limits for other charged sparticle from LEP are 
also imposed [56]. 

Constraints by dark matter relic density is also considered. In our analysis, the lightest 
neutralino is required to be LSP and a candidate of dark matter. The dark matter relic 
density is reported by WMAP7 in the range Qh"^ = 0.112 ± 0.0056 [57]. However, we only 
require the thermal abundance of neutralino satisfies a 3a upper-bound ^^oh^ ^ 0.1288 due 
to the reason that the neutralino might not be the only dark matter particle or that the 
neutralino might be produced in early universe via the so-called non-thermal mechanism 
[58]. The recent direct search experiment XENONIOO [59] is realized, which put the most 
stringent bound upon the dark matter-nucleon spin-independent scattering. We demand 
that the reduced dark matter-nucleon scattering cross section a^^ = a'^o 

^{n^ohynh^) should 

be smaller than XENONIOO limit. All constraints on flavor physics from low energy colliders 
and all bounds on dark matter from astrophysics and direct detection are implemented by 
using MicrOMEGAs [60] which uses the SuSpect output as input parameters. 

C. Numerical results and sparticle masses 

In this section we present some results based on a ~ 10^ case study in our parameter space 
scan. We find that most of the points can not yield correct RGE solutions nor induce the 
spontaneous electroweak symmetry breaking. Only ~ 700 points can pass all the constraints. 
To analyze the features of dark matter and collider signatures, it is found that light sparticles 
play the crucial roles. It is convenient to categorize the parameter points in term of the mass 
hierarchical relation of light sparticles as suggested in Ref. [40]. Here we define our four 
mass patterns: 

(1) the stop pattern (SO): m^o < m^^ < mf^,m^± ; 

(2) the stau/stop pattern (SS): m^o < m^^ < m^^ < m^± ; 

(3) the stau pattern (SA): m^o < m^^ < m^± < m^^ ; 



9 



(4) the chargino pattern (CH): m^o < m^± < mfj,mj^. 
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FIG. 1: Four mass patterns in the planes (a) m^o vs m^^ (top left), (b) m^± vs (top right) are 
illustrated. 



In Fig. 1, we present two scattering plots to reveal some features of these four mass 
patterns. In the plot (a), points are shown in the m^o — m^^ plane. In the plot (b), they 
are shown in m^± — m^^ plane. The distribution of these four mass patterns in soft SUSY 
breaking parameter space are also shown in Fig. 2. In Fig. 2(a), Fig. 2(b), Fig. 2(c), 
and Fig. 2(d), the distribution on the mo — M1/2 plane, on the mo,3 — Mi/2,3 plane, on the 
Aq — tan /3 plane, and on the Mi — fi plane are displayed, respectively. Four comments on 
the distributions are listed in order: 

(1) The gaugino mass running behavior in the non-universal scenarios is similar to that 
in mSUGRA. The 1-loop RGEs for gaugino masses are 



dt 



(4) 



where (&i, 62,^3) = (33/5,1,-3). It is well-known that in the mSUGRA the soft breaking 
parameters = M° = at Mqut evolves to M3 : M2 : Mi ~ 03]^° : OsM^ : aiM^ ~ 
6M3 : 2M2 : Ml at Mz- In the non-universal scenarios, gaugino mass parameters of first 
two generations at low energy hold this relation Mi : M2 ~ 1 : 2, which consequently means 
that the main component of the lightest neutralino x? is either Bino or Higgsino, depending 
on the relations of Mi and /i. This Mi — M2 relation also implies the lighter chargino Xi 
and the next lightest neutralino X2 should almost degenerate, since both of them are SU(2) 
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gaugino dominant m-+ ~ m^o ~ M2. In the Fig. 1(a) and Fig. 1(b), we can observe that 
the m^o and vary in the region of ~ (100, 350)GeV and (100, 700)GeV respectively, 
which basically reflects such a relation. On the other hand, the gluino mass parameter M3 
gains a large contribution via a negative (3 function at low energy, which leads to large gluino 
mass. It is beyond the reach of Tevatron, or even the reach of LHC with ^/s = 7TeV if 
gluino mass is heavier than 1 TeV. 

(2) The 1-loop RGEs for the third generation right-handed sfermion squared-mass param- 
eter could be written as 

16-^|< = %,^(m^„+mj3+m^3)-|,3^|M3r-|,,Wir + 4|Ar,r-l,?5, (5) 
16vr^|^l3 = %.^(<+<+^l3)-f^3W3p-^^7iWiP + 4|A,nP + ^^??5, (6) 
16vr2^m^-3 = 2y^,{m%^ + ml^ + ml^) - '^-^g\\M^\^ + '2\A,Yr\^ + ^-g\S, (7) 

where S is defined as 5* = Tr[yjm|.]. For the first two generation sfermions, the terms 
proportional to Yukawa couplings can be neglected due to the tiny values of their Yukawa 
couplings. In contrast, the sfermion mass terms of the third generation can get a large 
contribution from the Yukawa coupling terms. Therefore a light stop is often the lightest 
squark due to its large Yukawa couplings. When compared with mSUGRA, the assumption 
mo,3 < mo is can further decrease the contribution of gluino to the squark mass parameters. 
This can be read out from Fig. 2(b) where typically our stop pattern corresponds a smaller 
value in mo,3 parameter. The lightest squark is the lighter stop means that its cross section 
can be the largest at the hadronic colliders, Tevatron and LHC. The mass relation between 
stop and stau is difficult to be read out from RGE. Typically, the m? gets less positive 
contributions from gauginos and the Yukawa coupling is much smaller than Yt as well, 
which leads to quite a few mass patterns of light third generation sfermions. 

(3) The 1-loop RGE for the Higgs mass parameters if^ is given as 

leTT^^m^, = ^yl{rn\^ + + 772^3) - ^gl\M-,f - y^M^f + ^\at\^ + "^-gl . (8) 

The rnfj^ evolves to be negative at the low scale due to the large terms proportional to Y^^ 
which is essential in order to induce a spontaneous electroweak symmetry breaking. In the 
non-universal scenario, a smaller third generation sfermion mass leads to a |m|^^ | smaller 
than its counterpart in the mSUGRA. A small |m|^ | can lead to a small soft Higgs mass 
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FIG. 2: Four mass patterns in the planes (a) uiq vs M1/2 (top left), (b) mo,3 vs Mi/2,3 (top right), 
(c) Aq vs tan/3 (bottom left), (d) Afi vs fj, (bottom right) are displayed, respectively. 

parameter /x, due to the fact that /i is determined by the fnjj^ from the tree level relation 

„ Tnjj sin^(5 — Tnjj Cos'^lS mi , , 

^ ~ C0S2P 2 • 

When the magnitude of is comparable to that of Mi, as demonstrated in the cases of CH 

patten, the lightest neutralino Xi can have a large component of Higgsino due to the large 

mixing. 

(4) The squared-mass matrix for the stop quark in the weak interaction eigenstate basis 
{tL^in) is given by 

+ ml + A^^ v{alsin(3 — fiytCos(3) 

v{atsin/3 - fi*ytCos/3) m%^ + m^ + A^^ 
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where A<^. = (T^^ — Q^sin'^9w)cos(2l3)m%, and at = AtYf. The off-diagonal elements can 
generate a large mass split between two stop mass-eigenstates, which are labeled as (^1,^2)- 
For the case sign{(3) > chosen in this study, a small tgP and a large —at can produce a 
large mass splitting and leads to a relatively lighter stop ti, which is reffected in Fig. 2(c). 
In contrast, if the at is positive and the /i is small, m^-^ might be not light when compared 
with m^^, which leads to the cases of the CH patten. 

D. Features of the dark matter 

In SUSY models, most parameter space leads to a neutralino relic density which is too 
large to overdose the Universe. Depending on the mass spectra, four working processes can 
be introduced to produce a small relic density of neutralino [61]:(1) all the sfermions are light, 
neutralinos annihilate via t-channel sfermions exchange; (2) Xi has significant component 
of Higgsino or wino, main annihilating channel is to gauge bosons or Higgs; (3) neutralinos 
scatter with sfermions with nearly mass degeneracy which is so-called "co-annihilation"; (4) 
neutralinos annihilate via s-channel Higgs resonance with 2m^o = rriAo , or rriho , rriHo. 

When the stop ti is light, two main processes can lead to a small neutralino abundance. 
The first is the neutralino-stop coannihilation process Xi^i tg/ho, and the other is XiXi ~^ 
tt by the t-channel via exchanging the light stop when the kinematic is allowed. When the 
lighter stau fi is light, the neutralino-stau coannihilation can occur and make dominant 
contribution. This is the cases of SA and SS. It is worthy of mentioning that the neutralino- 
stop coannihilation can be significant for the case of SS. 

For the case of CH, the lightest neutralino Xi has a large component of Higgsino due 
to the large mixing deduced by the small fi and the large mo, as shown in the Fig.2. Two 
dominant annihilation processes can occur. (1) The pair of neutralinos can annihilation 
into gauge boson and Higgs boson via the processes XiXi W~^W~ , ZZ, Zho due to large 
Higgsino component. (2) A neutralino and a chargino can coannihilate into the particles of 
the SM via the processes XiXf ~^ W'^Z/'j, ff', which can occur due to the small chargino 
mass and a large Higgsino component in chargino xf- 

The last but not the least, the sfermion-sfermion self-annihilations processes may also 
be important, which can be attributed to the fact that a large sfermion-sfermion self- 
annihilation can increase the effective total dark matter cross section when the mass splitting 
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between sfermion and neutralino is small. Such a case can happen as we can read out from 
the effective cross section at the dark matter frozen-out epoch [62] 
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FIG. 3: The elastic spin independent DM-nucleon cross section cr^Op (left panel) and reduced cross 
section a^o ^ (right panel) for our four mass patterns are 



The neutralino-nucleon spin-independent scattering can be detected by dark matter di- 
rect search through XiQ ~^ XiQ via squark or Higgs exchange. If the first two generation 
squarks are light or neutralino has a large Wino/Higgsino component, the null results from 
direct detection can put a strong constraint on SUSY models. Recently, the XENON 100 
collaboration observed 3 events with an background of 1.8 ± 0.6 after 48kgx 100.9 days run- 
ning [59]. The result can be used to constrain the dark matter-nucleon spin-independent 
scattering cross section. 

We show the neutralino-nucleon cross section a^Op (left panel), the reduced cross section 
cr^Op (right panel) and XENONIOO limit in Fig.3. Here the reduced neutralino-nucleon cross 
section defined as c-Op ,, = a.^Op{Q^oh'^ /Qh"^) takes into account that the neutralino may only 
contribute part of the total dark matter relic density. From Fig.3 (a), we can observe that 
the neutralino in the CH case has a large a^Op due to its large Higgsino component. If 
we assume that the neutralinos are produced by a non-thermal process with a correct relic 
density Q^oh"^ = Qh"^, such a scenario is almost excluded by XENONIOO. If we assume the 
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neutralino only contributes part of the total dark matter relic density, such a scenario is still 
allowed, as shown in Fig. 3(b). 

IV. THE SIGNATURE OF LIGHT STOP PAIRS AT THE LHC 
A. production and decay of stop pairs 

In this section, we focus on the production and decays of the lighter top squark at the 
colliders. As pointed out above, the lighter stop ti can have a smaller mass than all other 
colored sparticles in the non-universal scenario and the mass splitting between it and the 
other colored sparticle can reach to 1 TeV. Then it is expected that the cross section of the 
lighter stop pairs at the LHC should be the largest one in our mass patterns introduced 
above. A pair of stop can be produced via two main processes via qq annihilation and gluon 
fusion [63]. The first one is dominant at Tevatron while the second is dominant at LHC. 

To evaluate the cross section of stop pair production a^^j^ at hadronic colliders, we utilize 
the package Prospino [64] which has incorporated the next leading order corrections. The 
results for Tevatron with ^/s =1.98 TeV, LHC with ^/s =7TeV and 14TeV are plotted in 
Fig.4. The K factor is determined to be situated in the range of ~(0.10, 0.13), ~(0.15, 0.20) 
and ~(0.15, 0.18) for Tevatron, LHC7 and LHC14, respectively. 

In the Fig.4, we observe that the mass of ti almost monotonically determines the cross 
section cij^, which is not sensitive to other SUSY parameters. At Tevatron, it is O.lpb or 
so for m^-^ ~ 200 GeV and quickly decreases to Ifb when m^^ increases to 400 GeV. At the 
LHC with ^/s =7TeV, the 0"^^^^ is larger than its value at Tevatron by a factor of 100 and 
decreases to Ifb for m^-^ ~ 800GeV. For the m^-^ <350GeV, more than one stop pair events 
at the LHC7 with 35pb~^ is predicted. With the increasing of ^/s to 14 TeV and a larger 
integrated luminosity (say 100 fb^^), stop pair events should be copiously produced when 
the lighter stop is less than 400 GeV. 

In order to analyze the signatures of stop pairs and to reduce the standard model back- 
ground by setting cuts, we have to study the decay products of stop pair. For this purpose, 
we use SDECAY package [65] to compute stop decay branching fractions. 

When the lighter stop is much heavier than the LSP, the two body decay ti — tx^ and 
ii bxt can be its dominant decay modes. In the ii bxt decay mode, the Xi can 
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FIG. 4: The cross sections of stop pair production at Tevatron, LHC-7TeV and LHC-14TeV are 
shown. 

mainly decay into Xi^~^ '^^ which is determined by the mass of stau and the dominant 
component of Xi ■ This is the feature of the SA and CH cases. When the hghter stop is 
smaller than 400GeV, these two body decay modes may be forbidden kinematically. Then 
the three body decay modes ti — )■ hWx\ and ti — )■ vfx^ become dominant ^ in our scan. 
When the mass splitting between stop and neutralino is too small to allow the three body 
decay, the four body decay mode ti — )■ bjjxi/bulxi (a recent study showed that this decay 
mode can be used to probe trilinear coupling Aq when tan/3 is small [66]) or the loop induced 
FCNC decay ii — )■ cx? become dominant. 

B. Simulation and analysis 

Recently, CMS and ATLAS collaborations reported several results on their searching for 
SUSY in different channels. In this section, we present the LHC bounds to the stop pair 
production with a collision energy ^/s = 7TeV in term of different channels studied by CMS 
and ATLAS. 

In our analysis, parton level events are generated by MadGraph [67], while parton shower, 
decays, and hadronization are performed by PYTHIA [68]. PGS [69] is used to simulate 

^ It is need to mention that the three body decay channels ii — > li^Xi searched by Tevatron coUaborations 
is not important in our scenario. Since the ruf is much smaher than mc,-, we observe that the branching 
fraction of ti vtx\ is always be larger than that of ti IC'Xi 
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detector effects and to find jets, leptons, and missing transverse momentum. The acceptance 
cuts for all jets and charge leptons are chosen as pt > 20GeV and |?7| < 2.5. 

The SUSY search strategies of CMS and ATLAS are optimized for mSUGRA scenario 
at this stage. It always require large MET and energetic leading jets in order to capture 
the signature from the heavy squark or gluino pair productions. Therefore the search to the 
jets plus MET channel sets the most stringent constraint to ordinary parameter space in 
mSUGRA. Cuts for studying this channel are described as follows: 

• In the Ref. [1], CMS presents a search for SUSY signatures on an integrated luminosity 
of 1.14fb"^ with jets and significant MET and without leptons. In this analysis, jets 
are required 

E^^''^' > lOOGeV, > 50GeV, \7]\j < 3, (12) 

where we have used the convention > The following selections are adopted 
to compute event rate 

Ht > 275GeV, ar > 0.55, (13) 
where Ht and are defined as 

Ht = T,E'^, ^T = \-^pi'\ (14) 

ar = E'j?I^H^-Ifl (15) 

If the jet number is more than 2, we utilize hemisphere algorithm [71] to combine jets 
into two pseudo-jets named Ji, J2 (assumed p;^^ > p;^^ again). The events containing 
isolated leptons and photons are rejected, where isolated leptons and photons are 
objects with 

p^T > lOGeV, pI > 25GeV, \ri\ < 2.5. (16) 

CMS collaboration performed a search in eight bins of Ht > 275GeV, and found the 
standard model background per bin could fit data well. The backgrounds decrease 
with increasing H-r rapidly, here we consider two bins 575 GeV< H^ <675GeV and 
Ht > 675GeV, in which the number of expected background is smaller than 20. We 
give the final event number in Fig. 5. 

• In the Ref. [2], ATLAS takes into account the events with jets and fJr with 1.04fb~^ 
of data. Events containing muons with px >10 GeV, |?7| < 2.4 and electrons with 
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FIG. 5: The event numbers of signal in the two bins given by the CMS jets + ^ searches with 
575 GeV< Ht <675GeV(left panel) and Ht > 675GeV(right panel) are shown. 



Pt >20 GeV, \rj\ < 2.47 are rejected. Reconstructed jets are required to satisfy 

p^^ > UOGeV, p'rj. > mCeV, \r]\j < 2.8. (17) 

The missing transverse energy ^ and the minimal frp) between frp and leading 

jets are required 

> ISOGeV, A0(j, prUin > 0.4 . (18) 

To probe the SUSY parameter points with different gluino and squark masses, the 
collaboration define five signal regions as 

A "2-jet" region, where cuts are chosen as nj > 2, me// >1000 GeV, ^/nieff >0.3; 

B "3-jet" region, where cuts are chosen as rij > 3, rrieff >1000 GeV, pr/meff >0.25; 

CI "4-jet" region, where cuts are chosen as rij > 4, mg// >500 GeV, ^/nieff >0.4; 

C2 "4-jet" region, where cuts are chosen as rij > 4, mg// >1000 GeV, ^/nieff >0.25; 

D "High mass" region, where cuts are chosen as Uj > 4, p^'-''" > SOGeV, rrieff >1100 
GeV, f^r/m^ff >0.2. 

Here rij is the number of the jets reconstructed in the event, the effective mass is the 
scalar sum of transverse momentum of jets and the transverse missing momentum 
The signal regions A, B, C and D correspond to different sparticle production channels 
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FIG. 6: The event number of signal after imposing the ATLAS lbjet+jets+f}p cuts: "2-jet" region 
(top left panel), "3-jet" region (top right panel), "4-jet high mgjj" region (bottom left panel), and 
"High mass" region (bottom right panel), are shown, respectively. 

qq, qg, light-^^ and heavy-^^, respectively. We show the final event number of signal 
passing all cuts in Fig. 6. The excluded values of sparticle production cross section 
are 22fb, 25fb, 429fb, 27fb, and 17fb, respectively. 

We also implement the bounds to the stop pair events in term of cuts from the SUSY 
searches at LHC with 35 pb^^ of data. The following channels have been included in our 
study : (l)CMS' jets + ^ channel [72], (2)ATLAS' jets + ^ channel [73], (3) ATLAS' 
llepton + jets + :^ channel [74], (4) ATLAS' 2leptons + jets + ^ channel [75], (5) ATLAS' 
lhjet + jets + IfjT channel [76], and (6) ATLAS' Ibjet + leptons + jets + ^ channel [76]. For 
the channels contained leptons, we do not focus on the lepton sign and flavor, and sum all 
the allowed events together. The bounds to the event numbers of the stop pair production 
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are shown in the Fig. 7. For these six channel, the upper hmits for new physics event 
number set by collaborations are 13.4, 45.5, 4.7, 20.7, 10.4 and 4.7 respectively. 

When comparing the bounds from 35 pb~^, we notice that the updated LHC bounds 
with 1 fb~^ of data could not further constrain on those non-universal models allowed by 
the bounds from 35 pb~^ of data. This is simply due to the fact that these cuts are not 
optimized to put bounds to the light stop pair production by the experimental collaborations. 
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FIG. 7: The event number of signal passed all cuts from LHC searches with 35pb ^ of data (from 
left to right and top to bottom): (l)CMS jets + fk, (2)ATLAS jets + (signal region A), 
(3)ATLAS llepto7i + jets + fr, (4)ATLAS 2leptons + jets + pr, (5) ATLAS Ihjet + jets + 
and (6) ATLAS Ibjet + leptons + jets + ^ are shown. 



C. Benchmark points 

Based on the mass pattern proposed above, we select four bench mark points (BMPs) 
to demonstrate their spectra and features at LHC. We tabulate the stop quark and LSP 
masses in Table I and provide the cross sections of these four bench mark points with the 
collision energy 7 TeV and 14 TeV. 

We focus on the dominant decay channel of the stop, which are listed below: 

BMPl : the stop dominantly decay via the t — tx^ channel {Br{t — )■ tx^) = 98.1%); 
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Benchmark points 


BMPl 


BMP2 


BMP3 


BMP4 


ruf-. 
m + 

Xi 

Al 


390 
207 
383 
206 


243 
471 
424 
223 


264 
199 
356 
190 


338 
179 
337 
176 


cj at 7 TeV (pb) 


0.23 


3.74 


2.33 


0.55 


a at 14 TeV (pb) 


2.54 


28.42 


18.91 


5.46 



TABLE I: The cross sections of pp — t- tit* in four bench mark points at LHC for colhsion energy 7 
TeV and 14 TeV are presented, respectively. The masses of hght sparticles (in GeV) are also given. 

BMP2 : the stop dominantly decay via the t cx° channel {Br{t cx°) = 98.7%); 

BMP3 : the stop dominantly decay via the t — )■ bu-rf hr + ^ channel {Br(t — j- brf]) = 
96.9%); 

BMP4 : the stop dominantly decay via the t -> 6 channel {Br{i bWx°) = 98.2%); 

Below we show several salient kinematic observables for these 4 bench mark points in 
Figs. (8-9). These observables offer us important clues as how signal can escape the current 
LHC SUSY search. The distributions shown in Fig. (8-9) explicitly show that the current 
cuts of ATLAS and CMS which are optimized for the mSUGRA search can not separate 
signal events of our bench mark points from the SM background events. There are several 
comments in order: 

• The reconstructed transverse missing momentum of these bench mark points is max- 
imal in the region with ^ < lOOGeV, as demonstrated in Fig. (8). Therefore the 
large missing energy cut adopted in experimental collaborations can greatly reduce 
the signals. 

• The effective mass for the bench mark point 2 and 4 is maximal in the region with 
m^ff < SOOGeV, while for bench mark point 1 and 4 the frieff can be large enough 
due to the large stop mass, as shown in Fig. (8). 

• The observables and mT^ are shown in Fig. (8). The observable is supposed to 
suppress QCD background heavily. Similarly, the razor method [77] , which is designed 
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FIG. 8: The transverse missing momentum, the effective mass of the four bench mark points, the 
aT and the mT2 quantity as well are demonstrated. 



to suppress the huge QCD background at LHC environment and to pick out the signal 
events from heavy pair-produced particles decay, can not help to distinguish signals 
from these bench mark points. While the reconstructed mT^ observable [70] is not 
large for all these bench mark points, the survival rate of signal after the cut with 
mT^ > 300 GeV can not be larger than 10% for bench mark point 1 and 1% for bench 
mark point 2, respectively. 

• The ratio of missing energy over Ht is shown in Fig. (9). The experimental cut on 
this quantity > 0.3 can affect bench mark point 1 and 4 significantly. 

• The transverse momentum of the leading two jets for the bench mark point is maximal 
in the region with Pt < 50 GeV region for bench mark point 2 and 3, as demonstrated 
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FIG. 9: The ratio of missing energy over the Ht, the angle separation of leading two jets AR{ji, j2), 
the transverse momentum of the leading jet, and the number of jets distribution in each bench 
mark point are demonstrated. 



in Fig. (9). Therefore the cut demanding the leading jet must be larger than 150 
can affect signal significantly. Meanwhile, such a cut can also considerably reduce the 
signal from bench mark point 1 and 4. Similarly, the cut on the second leading jet can 
affect all bench mark points badly. 

Next, we focus on the signature from bench mark point 1 and 2, and consider the corre- 
sponding SM backgrounds. The final states of these two bench mark points are provided in 
Table (II) and Table (III), where we have imposed the acceptance cuts to both of them: 1) 
Ptij) > 20 GeV, 2) Pt{i) > 20 GeV, 3) ^> 20 GeV, 4) r/(j) < 2.5, and 5) R{i,j) > 0.4 and 
> 0.4 as well. 

Let's first look at the signal from the first bench mark point. The search strategies for 
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2j 


3j 


4j 


5j 


>6j 


ne = 


3% 


8% 


15% 


16% 


19% 


ni = l 


3% 


7% 


7% 


4% 


2% 


n£ = 2 


0.9% 


0.6% 









TABLE II: The percentage of lepton and jet multiplicity channels determined by our bench mark 
point 1 are shown. 





2j 


3j 


4i 


5j 


> 6i 


ne = 


20% 


9.5% 


3.0% 


1% 





TABLE III: The percentage jet multiplicity channels from the bench mark point 2. 

searching ti tx^ can be categorized by the final states in the literature: the full hadronic 
channel [78, 79], the semi-leptonic channel[80, 81], and the di-leptonic channel. For the 
hadronic and semi-leptonic modes, some kinematic observables have been studied to separate 
the signal and background, i.e. the missing energy and effective mass. Moreover, When the 
tops in the final state are highly boosted, the top tagger based on the jet substructure 
analysis can be used to distinguish signal and background [82, 83]. It is also remarkable 
that due to the right-handed helicity of the light stop, the top quarks in the final state 
should be polarized [84]. The di-leptonic channel pp — tit^ tix^X^ ~^ ^^bb + is less 
studied in literature, the apparent reason is the small branching fraction. However, consider 
the messy background at LHC, two leptons in the final state can help to suppress background 
greatly. Furthermore, in order to claim the signal is from light stop pair decay, di-leptonic 
channel should also be observed. So di-leptonic channel is complementary for the discover 
of light stop and deserves careful study. 

Now, let's look at the signature of the second bench mark point. Search for the signature 
ti — > cxi at LHC also existed in a vast of literature. If ti dominantly decays into charm 
and neutralino, it might be quite challenging to directly detect it even it might be copiously 
produced, similar decay could occur in the T-parity little Higgs model as demonstrated 
in [85]. The charm jet can be very soft which might escape the triggering, meanwhile 
the reconstructed missing energy can not be large. Therefore such events might not even 
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be recorded. It is known that b-tagging and top-tagging can increase the supersymmetry 
signal relative the standard model backgrounds. Therefore, c-tagging techniques are also 
suggested to distinguish signal and QCD background [86, 87]. 



Distribution of m^,, with signai and background stacked The distribution of with signal and bacltground stacked 




GeV GeV 



FIG. 10: The distributions of effective mass rueff and the reconstructed transverse mass for 
the channel f]+ i + jets after ah other ATLAS cuts for bench mark point 1 in 7 TeV are shown. 




In Fig. (10-11) we further show how the current ATLAS cuts and the current CMS cuts 
affect the observation of our bench mark point 1 and 2, respectively. The QCD background 
in Fig. (11) can be safely omitted due to the fact that both the ax and Rmis can suppress 
it significantly. The plots show the distribution of the kinematic observables rrieff and rri^ 
with the current ATLAS cuts and m|. and me// with the current CMS cuts. After imposing 
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all cuts, we arrive at the results given in Table IV. The results indicate that these bench 
mark points, especially for bench mark point 2, define a challenge to the current searching 
strategy. For bench mark point 1, we have studied the full hadronic channel and the di- 
leptonic channel (mainly suppressed by the branching fraction) and the bound from the 
semi-leptonic channel is more stringent. 

Although the cross section of signal increased by a factor 10 or so from 7 TeV to 14 
TeV, the cross section of background also increase with almost the same factor. We arrive 
at the conclusion that even when LHC can run with the collision energy 14 TeV, it is still 
challenging to explore these bench mark points. 





signal 


background 


S/B 


S/^/S + B 


Lum. (7 TeV) 


BMPl 


0.04 


4.0 [74] 


1 X 10^2 


0.02 


62.5 fb-i 


BMP2 


0.01 


24.5 [72] 


4 X 10"^ 


0.002 


6250 fb-i 



TABLE IV: Number of events after CMS and ATLAS search cuts with integrated luminosity 35 
pb~^ are displayed and the required luminosities with the collision energy 7 TeV for the discovery 
S/ \/ S + B = 5 of bench mark points 1 and 2 are estimated. 

Signature of bench mark point 3 and 4 is bbT~^T~, where both b jets and the rs are 
typically soft due to the small mass splittings in the cascade decay chains. A recent study 
on such final states in mSUGRA context can be found at [88]. The b-tagging and tau-tagging 
performance with a soft b jet and a soft tau decreases which may challenge the success to 
distinguish signal events. Similar to the signature of bench mark point 1 and 2, the current 
cuts on rrif^ff and -^suppress the signal badly, therefore to find signature of these two bench 
mark points at LHC is also difficult. We neglect the detailed analysis of them. 

D. Search at the ILC 

The ILC is a future electron-positron collider. At the first stage, it will start with 220 
GeV and run up to the maximal center of mass energy 500 GeV. At its later stages, it can 
be upgraded up to 1 and 3 TeV. In principle, it is also designed to be able to scan near the 
threshold region of particle production. Compared with hadronic collider, the ILC enjoys 
much clean environment as well as high energy resolution capability. For the discovery of 



26 



SUSY, as one of its advantage, it can reconstruct three dimensional momentum of missing 
energy. 

From the ILC detectors [89-91] design report, we take the energy resolution for jets is 
assumed as 



Meanwhile, since the detectors are supposed to cover the full solid angle, we assume that 
the coverage of detectors to charged tracks can reach to 20 for rj in our fast simulation. Such 
a detector simulation is realized by modifying the PGS card file. 

At the ILC, a pair of light stops can be produced via the s-channel 7/Z exchange. As 
pointed out in the design report [92], ILC can cover the region of parameter space with the 
light stop. Several realistic Monte Carlo study have been performed in the literature [93-98]. 
Below we focus on the detection of the signal represented by our bench mark point 2 [13]. 

For the bench mark point 2, the pair production of neutralinos should be possible but 
can not overcome the background e'^e" — t- uu. Then the light stop should be the first super 
particle detected at ILC. Here we update the relevant analysis by considering more kinematic 
observables and using the neural network discriminant to improve signal and background 
separation. And we find that when compared with the sequential cut method, the neural 
network discriminant analysis can improve both the ratio of signal over background and the 
significance remarkably. Our analysis can be extended to the similar study for b — )■ bx^ and 
we expect the ratio of signal over background and the significance can also be considerably 
improved. In Figs. (12), we show the distribution of key kinematic variables used as the 
input of the neural network analysis, where we do not impose any a cut except the transverse 
momentum of jets Ptij) > 5GeV: 

• the reconstructed missing energy fJ, which is obtained from ^/s — Ej^ — Ej^. This 
quantity can not be reconstructed at LHC but can be reconstructed at ILC. For signal 
this quantity should be large, as shown in Fig. (12). While for the background ZZ 
which is occurred in t and s channel, when one of Z decays invisibly, the missing energy 
can be around 250 Gev, this explain the bump in the ZZ background. 

• The ratio of ^/rrieff, where Tn^ff is defined as the visible energy sum of all objects 
in the event. Obviously, for signal, this quantity should be large, as shown in Fig. 
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FIG. 12: The distributions of the reconstructed missing energy, the invariant mass of leading two 
jets, the ratio of ^/ Ej^ , and the invariant mass of missing Lorentz vector as well are demonstrated. 
The unit of the y-axis is determined by normalizing the integrated luminosity to lfb~^. 



(12). From the distribution of this quantity, it becomes quite clear that the dominant 
background events are from WW and ZZ . 

• The centrality C. In the signal, the energy tends to deposit in the direction with 77 = 
region. We find this quantity is useful. 

• The transverse momentum scalar sum Ht and the jet mass of the two hemisphere jets. 
If there are more than 2 jets, we can use the hemisphere algorithm to group jets into 
two fat jets. For signal, the invariant mass of each jet should be small. While for 
background, like the highly boosted weak bosons, the invariant mass must be large. 
We find that these observables are useful to suppress background. 
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We can choose a few pre-selection rules and adopt a few simple cut method to suppress 
background while maintain a good acceptance to signal. At the pre-selection level, we use 
lepton veto and a cut on ^> SOOGeV". Then the dominant background after pre-selection 
is WW pair and eeZ events. We also list the results in the simple cut method, where we 
choose: 1) > 425GeV, 2) p/{Ej^ + Ej^) > 10, and 3) m(ji,j2) < 60. The results are 
presented at the third line of Table V. Due to the correlation among kinematic observables, 
it is difficult to find the best set of cuts. To finish such a task, we utilize the neural network 
discriminant analysis to optimize cuts. The results are presented in Table (V). 





signal 


tt 


WW 


euW 


zz 


eeZ 


S/B 


S/VS + B 


No. of Events after pre-selection 


11.1 


6.2 


336.7 


8.9 


44.8 




0.03 


0.54 


No. of Events after a few cuts 


11.1 




18.6 


1.0 


0.7 




0.5 


1.9 


No. of Events after NN 


9.6 




0.9 


0.7 


0.4 




4.8 


2.6 



TABLE V: The number of events ( normalized to the integrated luminosity 10/6 ^ ) after prese- 
lection, after some simple cuts and after the NN discriminant cut are demonstrated. 

The results of the neural network discriminant are presented at the forth line of Table V. 
It is obviously that the NN discriminant analysis can improve both the ratio of signal over 
background and the significance. From the results of the neural network analysis, we can 
estimate the required luminosity is 37 fb~^ for the discovery significance S/yJ S + B = 5. 

A simpler version of the neural network discriminant analysis can be demonstrated in 
Fig. (13), where two pairs of two almost independent observables are shown. The neural 
network discriminant analysis basically utilizes such the correlations among observables to 
distinguish signal and background. 

V. CONCLUSION AND DISCUSSIONS 

We have studied the non-universal SUSY models and explored light stop pair production 
at the LHC We scan the SUSY parameter space at the GUT scale and evolve to low energy 
considering the bounds by LHC searches with 35 pb~^ and 1 fb~^ of data and the dark 
matter relic density and direct search bounds. We find that to give correct relic density the 
stop usually has small mass difference with neutralino. Such a scenario easily escape the 



29 




FIG. 13: Two dimensional scatter plots for signal and background are demonstrated. The blue and 
red points denote signal and background respectively. In the left panel, the correlations between 
the rate and centrality for signal and background is shown. In the right panel the correlation 
between the invariant mass of the missing 4-momentum and the invariant mass of the two leading 
jets is shown. 

current search cuts adopted at both CMS and ATLAS collaborations. 

The model we explored in the work is an important scenario since the colored SUSY 
particles are the primary goal to search at LHC and gluino is usually very heavy when 
evolving from GUT scale to the low energy scale. Further we have to consider the dark 
matter relic density bound, at least considering the upper bound so as not to overdose the 
Universe. The dark matter relic density usually leads to degenerate pattern between the 
light stop and the LSP neutralino. 

In this work we demonstrate that it is difficult to detect the light stop scenario if only the 
pp — )■ titl process is considered. There have been some studies to improve the ratio between 
signal and background using the associate production, such as via associate mono-jet and 
mono-photon processes [16, 99]. It is found by using the associated production, the large 
region of the parameter space can be covered by LHC. Another method proposed to further 
suppress the background is by utilizing the two energetic tagged b jets [18]. By studying the 
pp — )■ iiilbb process [18], even the very degenerate stop neutralino scenario can be explored 
for LHC at 14 TeV and 500 fb~^ integrated luminosity. 

Although the associate production channels offer a hope to detect such a difficult scenario 
at LHC, the required luminosity seems too large. Our simulation shows that the ILC is an 
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ideal place to probe these models. Whether the SUSY is hidden at LHC or SUSY does not 
exist at low energy may need more careful and fortitude studies. 
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